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ABSTRACT 

We study the tidal disruption of binaries by a massive point mass (e.g. the 
black hole at the Galactic center), and we discuss how the ejection and capture 
preference between unequal-mass binary members depends on which orbit they 
approach the massive object. We show that the restricted three-body approx- 
imation provides a simple and clear description of the dynamics. The orbit of 
a binary with mass m around a massive object M should be almost parabolic 
with an eccentricity |1 — e| £ (m/M) 1 / 3 <C 1 for a member to be captured, while 
■ the other is ejected. Indeed, the energy change of the members obtained for a 

parabolic orbit can be used to describe non-parabolic cases. If a binary has an 
encounter velocity much larger than (M/m) 1 / 3 times the binary rotation velocity, 
it would be abruptly disrupted, and the energy change at the encounter can be 
evaluated in a simple disruption model. We evaluate the probability distribu- 
tions for the ejection and capture of circular binary members and for the final 
energies. In principle, for any hyperbolic (elliptic) orbit, the heavier member 
has more chance to be ejected (captured), because it carries a larger fraction of 
the orbital energy. However, if the orbital energy is close to zero, the difference 
between the two members becomes small, and there is practically no ejection 
and capture preference. The preference becomes significant when the orbital en- 
ergy is comparable to the typical energy change at the encounter. We discuss its 
implications to hypervelocity stars and irregular satellites around giant planets. 

Subject headings: binaries: general, Galaxy: kinematics and dynamics, planets 
and satellites: formation, Galaxy: Center, Galaxy: halo, Planets and satellites: 
individual (Triton) 
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1. Introduction 

The disruption of a star by a massive black hole (BH) is one of the most spectacular 
examples of the tidal phenomena (Komossa & Bade 1999; Donley et al. 2002; Grupe et al. 
1995). A star that wanders too close to a massive BH is torn apart by gravitational forces. 
Almost half the debris would escape on hyperbolic orbits, while the other half would traverse 
elliptic orbits and returns to periapsis before producing a conspicuous flare (e.g. Rees 1988). 
The disruption process has been numerically investigated in detail (Evans & Kochanek 1989; 
Laguna et al. 1993; Ayal et al. 2000; Kobayashi et al. 2004; Guillochon et al. 2009), and 
the new generation of all-sky surveys are expected to detect many tidal flares (Strubbe & 
Quataert 2009, Lodat & Rossi 2011). Recently a possible discovery of the onset of the rapid 
BH accretion has been reported (Burrows et al. 2011; Zauderer et al. 2011; Levan et al. 
2011; Bloom et al. 2011). 

Once a star gets deeply inside the tidal radius of a BH, the tidal force dominates over the 
self-gravity and thermal pressure of the star. A very simplified description of the disruption 
process could be the encounter between a star cluster (or a cluster of point masses) and a 
massive BH. The simplest case consists of a binary and a massive BH in which after the 
tidal disruption, one star would escape to the infinity, while the other could be captured 
by the BH. This is actually one of the leading models for the formation of hypervelocity 
stars (Hills 1988; Yu& Tremaine 2003). The captured stars may explain the S-stars in the 
Galactic center (Gould & Quillen 2003; Ginsburg & Loeb 2006; Ghez et al. 2005; Genzel et 
al. 2010). 

Hypervelocity stars are stars with a high velocity exceeding the escape velocity of the 
Galaxy. After the discovery of such stars in a survey of blue stars within the Galactic halo 
(Brown et al. 2005; Hirsch et al. 2005; Edelmann et al. 2005), many authors have predicted 
the properties of the hypervelocity stars (Gualandris et al. 2005; Bromley et al. 2006; Sesana 
et al. 2007; Peretz et al. 2007; Kenyon et al. 2008; Tutukov & Fedrova 2009; Antonini et 
al. 2010; Zhang et al. 2011). These investigations so far have used three-body simulations 
or analytic methods that relied on results from three-body simulations. 

The six orders of magnitude mass ratio between the Galactic center BH and the binary 
stars allows us to formulate the problem in the restricted three-body approximation. In a 
previous paper (Sari et al. 2010, hereafter SKR), we have shown that the approximation 
is efficient and useful to understand how binary stars behave at the tidal breakup when 
the binary's center of mass approaches the BH in a parabolic orbit. In this paper, we 
generalize the approximation for orbits with arbitrary eccentricity. This enables us to give 
a complete picture of the ejection and capture process. We also provide the ejection and 
capture probability distributions that can be simply rescaled in terms of binary masses, their 



-3 - 



initial separation and the binary-to-black hole mass ratio when applied to a specific system. 
Our method is computationally more efficient than full three-body simulations, and it is 
easier to grasp the nature of the tidal interaction. 

In §2, we outline the restricted three-body approximation. In §3, we evaluate how much 
energy each member gains or loses at the tidal encounter and we discuss how the energy 
change evaluated for a parabolic orbit can be used to study non-parabolic orbit cases, and 
in §4, we give qualitative discussion on the ejection and capture preferences. In §5, we study 
high velocity encounters. In §6, the numerical results are discussed. In §7, we use our results 
to describe the capture process of Triton around Neptune. Finally, in §8, we summarize the 
results. 



2. The Restricted Three-Body Problem 

The equation of motion for each of the binary members is given by 

GM Gm 2 
r i = — ri + 1 ^(rs-n), (1) 

f2 = _£M r2 _^- (r2 _ ri) , (2) 

?2 |ri — r 2 \ 

where V\ and r 2 are the respective distance from the massive point mass with M. We will 
call the point mass the black hole (BH), though the binary is assumed to travel well outside 
the event horizon and our results can be applied to any systems which include a Newtonian 
massive point mass. The equation for the distance between the two r = r 2 — r\ is 

GM GM Gm 

-r 2 H q-ri r r. (3) 



r'2 r'{ r 3 

where m = mi + m 2 <C M. We assume that the two masses are much closer to each other, 
and to the trajectory of the center of mass of the binary r m , than each of them to the BH. 
Both energy and orbit obtained under the approximation are fairly accurate except a part 
of the orbit just around the periapsis passage (see SKR for the details). 

Linearizing the first two terms of equation fl3]) around the center of mass orbit r m , we 
find that the zero orders cancel out. Then, rescaling the distance between the bodies by 
(m/M) 1 / 3 r p and the time by ^Jr 3 / GM where r p is the distance of the closest approach 
between the center of mass of the binary and the BH, we can re-write eq. ([3]) in terms of the 
dimensionless variables: r] = (M/m) 1 / 3 (r/r p ) and t: 

f,= f^) 3 h»7 + 3(r;f m )r m ]-^, (4) 

\ ' m / I'll 
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where r m is a unit vector pointing the center of mass of the binary. We define the orbit of the 
center of mass to be a conic orbit r m / r p = (l + e)/(l + e cos /) where e is the eccentricity, and 
the true anomaly / is the angle from the point of closest approach. Since r m = (cos /, sin /, 0), 
and we set r) = (x,y,z), explicit equations in terms of dimensionless Cartesian coordinates 
reads 

(1 + ecos/) 3 . x 

x = — n i \3 — + 3 {xcosf + ysmf cos/ - — , , , 5 

(1 + ecos/) 3 . . y 

f = (1 + e )3 + 3 ^ C ° S f + y 5111 /) 5111 /l ~ ( X 2 + y 2 + ^2)3/2 » ( 6 ) 

(1 + ecos/) 3 g 

(1 + e) 3 (x 2 + y 2 + 2 2 ) 3 / 2 ' K ) 

where the eccentricity 

is related to the energy of the center of mass which is given by 

^ m.. , 9 GMm , , 

£ = 2-W 2 --— ■ (9) 

z, / m 

Using the dimensionless time, the conservation of the angular momentum can be expressed 
as 

/ = (l + e )" 3 / 2 (1 + e cos/) 2 . (10) 

Analytically one has relations between r m and t through a parameter which are given by 
(e.g. Landau & Lifshitz 1976) 

E<0 r m /r p = (l-e) _1 (l -e cos £), i = (1 - e)~ 3/2 (f - esin£) , (11) 
E = r m /r p = (1+e 2 ), t = V2(Z + £ 3 /3), (12) 
£>0 r m /r p = (e - (e cosh £ - 1) , t = (e - 1)" 3/2 (esinh^ - , (13) 

where the closest approach r m = r p happens at t — 0. 



3. Energy change at the BH encounter 

We are interested in the fate of stars in a binary, following its encounter with a massive 
BH. In order to study the ejection and capture process, we evaluate the energies of the stars 
as functions of time. When the binary is at a large distance from the BH, the binary members 
rotate around their center of mass which gradually accelerates towards the BH. The specific 
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self-gravity energy of the binary is about — Vq = —Gm/a. Analytic arguments (SKR) suggest 
that at the tidal breakup one member gets additional energy of the order of v m v where v m is 
the velocity of the center of mass at the tidal radius r t = [Mj m) 1 / 3 a. If the binary approaches 
the BH with negligible orbital energy, the velocity is v m = [GM /r t ) l l 2 = Vq(M /m) 1 / 3 . 
The additional energy is larger than the self-gravity energy by a factor of (M/m) 1 / 3 3> 1. 
Therefore, we will neglect the self-gravity term in the following energy estimates. This 
treatment is valid as long as the binary is injected into the orbit r m at a radius much larger 
than the tidal radius. 

The energy of one binary member is given by 

e, = !* N . _ 9Msa.. ( i4) 

Linearizing the kinetic and potential energy terms around the orbit of the center of mass r m 
and using the initial energy Jj = (rrii/m)E, we obtain 

Ei = Ii + AEu (15) 

. GMrrii 
AEi = mjr m (ri-r m )H - — r m (r. ; - r m ) , (16) 

Since in our limit the total energy of the system is E, considering AE 2 = —AEi, we get 

A E 2 = ( r mi . + ^r m r ) . (17) 

Using our rescaled variables, the additional energy is given by 

Gm x m 2 fM\ 1/3 - 

AE 2 = -AE X = — — AE, 18 

a \ m J 



AE = D- 1 



-xsmf + y(e + cosf) (1 + ecos/) 2 . ' 

H — [X cosf+y sin /) 



VTT^ (l + e) 



(19) 



where D = r v j r t is the penetration factor which is useful to characterize the tidal encounter. 
Once the binary dissolves, AE becomes a constant because the body is eventually moving 
only under the conservative force of the BH. Hereafter, the energy change Ai? means the 
constant value after the disruption, otherwise we specify it. The equation of motion (j3J 
indicates that the negative of a solution r = r(£) is also a solution. The energies AEi are 
also linear in the coordinates. Therefore, another binary starting with a phase difference 
7r will have the same additional energy in absolute value but opposite in sign. A uniform 
distribution in the binary phase implies that, when the binary is disrupted, each body has a 
50% chance of gaining energy (and a 50% chance of losing energy). 
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As we have discussed, the typical energy change is larger than the self-gravity energy 
by a factor of ~ (M/m) 1 / 3 , it is of order of (Gmirri2/a)(M /m) 1 / 3 . Then, the dimensionless 
quantity AE is an order-of-unity constant after the disruption. Its exact value depends on 
orbital parameters, but for qualitative discussion we just need to know that AE is about 
unity. Later we will numerically show that AE is an order-of-unity in the relevant parameter 
regime^, and numerical values will be used to estimate the ejection and capture probabilities. 

Rescaling energies by the typical value of the energy change, the energies of the binary 
members after the disruption are given by 

E x = h- AE, E 2 = I 2 + AE, (20) 

where bar denotes energy scaled by (Gmim2/a)(M /m) 1 ^ 3 . An interesting outcome of the 
encounter between a binary system and a massive BH is the "three-body exchange reaction" 
(Heggie 1975; Hills 1975) where one member of the binary is expelled and its place is taken 
by the BH, i.e. one binary member is captured by the BH and the other is ejected to infinity. 
In order for a member rrii to escape from the BH, the initial binding energy should be smaller 
than the energy gain: < \AE\ ~ 1. The same condition is required when a member of 
the binary in a hyperbolic orbit loses energy and it is bound around the BH. Therefore, when 
we discuss the ejection or capture process associated with a massive BH, the absolute value 
of the initial energy should be comparable or less than unit: \I^\ £ 1. 



Since the energy, penetration factor (periapsis radius) and eccentricity are related by eq 
(JEJ) or equivalently: 

only two of them are the independent parameters to describe the binary orbit. Considering 
\Ii\ 1 together with the mass ratio (m/M) 1 / 3 <C 1 and the tidal disruption condition D £ 1, 
the eccentricity should be almost unity |1 — e\ D{m/M) l / 3 {m par /m) for a member rrii to 
be ejected or captured where m par is the mass of the partner {m paT = m 2 for i = 1 and mi 
for i = 2). If we use the semi-major axis r a = r t D/(l — e), the condition can be rewritten 
as \r a /r t \ ~ (M/m) 1//3 (m/m par ) where r a is negative for hyperbolic orbits. 

Such orbits differ very little from parabolic orbits with the same periapsis distance, 
especially around the tidal radius and inside it. Therefore, the energy change AE is expected 
to be almost identical to that for the parabolic case. As long as we study the exchange 



1 For prograde orbits with D ~ 0.1, the energy change is as large as AE ~ 30 in a very narrow range 
of the binary phase (see figure 7 in SKR) where the binary members once come close to each other before 
they break up. However, the phase-averaged value (lAi?!) is still an order-of-unity and it is a more relevant 
quantity for the discussion on the ejection and capture probabilities and the final energies 
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reaction, we can approximate AE by the parabolic results AE e= i. However, e ~ 1 does 
not necessarily mean \E\ <C 1. In general, we need to take into account the offset of the 
final energy due to the non-zero initial energy, which would affect the ejection and capture 
probabilities. The final energies are approximately given by 

E 1 = I 1 -AE e=1 , E 2 = I 2 + AE e=1 . (22) 

4. Which gets kicked out? 

We here consider a simple question: Which member is ejected or captured if an unequal- 
mass binary is tidally disrupted by a massive BH? If E > (hyperbolic orbits), one binary 
member could lose energy and get captured by the BH, the other flies away with a larger 
energy. Assuming a uniform distribution in the binary phase, each member has a 50% chance 
of losing energy (and gaining energy). However, since the lighter one (the secondary) has 
a smaller initial energy, it is preferentially captured and the heavier one (the primary) has 
more chance to be ejected. 

For elliptical orbits, by considering a plausible semi-major axis r a , we can obtain tighter 
constraints on the eccentricity and energy, compared to the requirements from the exchange 
reaction. This is particularly relevant for studies of hypervelocity stars. If r a is around the 
radius of influence of the BH ~ GM/a 2 where a is the local stellar velocity dispersion, 
for the Galactic center, it is about r a ~ a few parsecs ~ 10 5 r t for a ~ several solar radii. 
Then, we get 1 — e = D{r t /r a ) ^ 10~ 5 and \E\ ~ [r t /r a ){M /m) 1 ^(m/mi){m/rri2) ~ 10~ 3 . 
Our previous estimates based on parabolic orbits are appropriate to study the production of 
hypervelocity stars for which an equal ejection chance is expected (SKR). When the semi- 
major axis is as small as r a ~ (M / m) 1 ^ (m / m par )r t , the initial energy would be of the 
order of unity as we have discussed and it affects the ejection preference. Since the secondary 
has less negative initial energy, it is preferentially ejected. 

Recently, Antonini et al. (2011) performed N-body simulations of unequal mass binaries 
with mi = 6M Q , m 2 = 1 or 3M and a = 0.1 AU in elliptical orbits around a supermassive 
BH M = 4 x 1O 6 M . They find that the initial distance of the binary from the central BH 
plays a fundamental role in determining which member is ejected: for a large initial distance 
d = O.lpc or equivalently r a ~ 3 x 10 3 r i; the ejection probability is almost independent on the 
stellar mass, while for d = O.Olpc or r a ~ 3 x 10 2 r t , the lighter star is preferentially ejected. 
Considering that the ejection probability significantly decreases if r a becomes smaller than 
~ (M /my^^/rripa^rt ~ 80(m/m par )r t , these results are consistent with our analysis. 

These ejection preferences for hyperbolic and elliptic orbits are naturally understood if 
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we consider a large mass ratio for the binary members. The energy of the primary practically 
does not change at the tidal encounter. Whether it is ejected or captured after the tidal 
breakup simply depends on the initial energy ~ E, while the secondary might have a chance 
to make a transition between a bound and unbound orbits around the BH (Bromley et al. 
2006). In the large mass ratio limit, the exchange reaction condition (i.e. the transition 
condition for the secondary) is \E\ ^ (Gm 2 / a)(M /m) 1 ^ 3 or equivalently |r | ^ (M/m) 2 / 3 a. 



5. High Energy Regime 

If a binary has a large orbital energy E ^> 1, both members are ejected after the BH 
encounter as a binary system or two independent objects 0. Although the high energy regime 
is not important in the context of the three-body exchange reaction, we discuss the regime to 
clarify the parameter dependence of the numerical results in the next section. A high orbital 
energy E 3> [M/m) l ^{m/mi){m/rri2) affects the velocity of their center of mass at the 
encounter v m ~ (E/m + GM/r m ) 1 / 2 = (M/m)^ 3 v y/(e - I)/ 2D + r t /r m . Then, the tidal 
disruption radius (i.e. where a binary is disrupted) can be defined in three different ways. 
We here order them from a large to small radius, (a) Relative acceleration: the radius at 
which the BH tidal force becomes comparable to the mutual gravity of the binary. This is r*. 
(b) Relative velocity: the radius at which the tidal force induces the relative velocity between 
the binary members comparable to the binary escape velocity Vq. (c) Relative position: the 
radius at which the difference in position increases by more than the initial binary separation. 

The duration that the center of mass is around r m is of order of At ~ r m /v m . During 
this period, the tidal acceleration of the relative motion of the binary members by the BH is 
of order of A ~ GMajr 3 m . The two radii (b) and (c) can be estimated from two conditions: 
Av = AAt ~ Vq and Ax = AAt 2 ~ a, provided that the duration of the encounter is 
comparable or shorter than the binary rotation time-scale: At ^ a/vo- If the energy is high 
(e — 1)/D = 2E(m/M) l l 3 (mi/m)(m2/rn) 3> 1, these conditions give r m = r t D 1 ^ 4 /(e — l) 1 / 4 
and r t D/(e — 1), respectively. Since they should be larger than the periapsis distance, only 
the cases that satisfy D ^ (e — l) -1 / 3 for the radius (b) or e ^ 2 for the radius (c) lead to 
the disruption. The radius (c) is basically the place at which the orbit of the center of mass 
makes its turn (i.e. the periapsis). If the energy is low (e-l)/D<l, all the estimates give 
the original tidal radius r t . 

When we discuss the energy change AE at the tidal encounter, there are two important 



2 If E is a large negative value, both members are captured after the disruption. Since the velocity at 
the tidal radius is reduced v m ~ (M/m) 1 / 3 -^, the energy change should be smaller \AE\ 1. 
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points which we should emphasize. First, the energy of each of the binary members in the 
BH frame changes only due to the mutual force between the binary members. Secondly 
most of the work done by one member on the other, which is AE, is done outside the 
tidal radius r t . The mutual force is of order of Gm 1 m 2 /a 2 . During the binary rotation 
time-scale a/vo, the force acts over a length ~ (v m /vo)a in the BH frame. Therefore, the 
work is W ~ (mim2/m)v m VQ. Since the direction of the mutual force changes with the 
binary rotation, AE(t) oscillates with the amplitude of W. When the binary is disrupted, 
AE becomes a constant value which is basically determined by the binary phase at the 
disruption. Then, we might expect that the final value of AE is a sinusoidal function of the 
binary phase for circular binaries. As we will see later, this is actually the case for the high 
energy encounters. Even with the largest estimate of the tidal radius (i.e. r t ), the duration 
of the encounter At ~ r t /v m is shorter than the binary rotation time-scale by a factor of 
~ sjD j[e — 1) <C 1. The work during the encounter is negligible compared to the work W 
which has been done outside r t . On the other hand, in the low energy regime, the duration 
of the encounter is comparable to the binary rotation time-scale. Considering that at the 
encounter the orbits of the members in the comoving frame of their center of mass should be 
significantly deformed from the original orbits (e.g. circular orbits) before they finally break 
up, the work during the encounter could induce deviation of AE((p) from a simple sinusoidal 
function. However, the typical value is still expected to be about AE ~ [m,\mij / m)voV m . 

In both the low and high energy regime, the typical energy change is given in a dimen- 
sionless form by 

where we have assumed r m = r t to estimate v m . In the high energy regime, the disruption 
might happen at a smaller radius, but v m is determined by the orbital energy and it is 
insensitive to the choice of r m . When (e — 1)/D 3> 1, the energy change becomes much 
larger than unity. However, the energy gain is not significant compared to the original 
energy E, and one finds that the tidal encounter is not an efficient acceleration process 
anymore. 

In the high energy regime, the energy change eq f fT9|) can be evaluated by assuming 
that a binary is abruptly disrupted at the tidal radius r t , since the work during the tidal 
encounter is negligible. For a circular coplanar binary: (x, y) = -D _1 (cos0t, sin0 t ) and 
(x,y) = ±D 1 / 2 (— sin t , cos t ), we obtain 

AE = ( 1 ± - 1 ) sin /t sin0 t +(l± 1 + 



COS ft COS ( 



(24) 
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where <\>t is the binary phase at the tidal radius, ft is the negative value solution of 1 + 
e cos ft = (1 + e)D and the signature indicates a prograde (+) or retrograde (-) orbit. This 
is a sinusoidal function of the binary phase as we expected, and the square of its amplitude 
is 3 — r t /r a ± 2-^/(1 + e)D which is larger for prograde orbits and the difference between 
prograde and retrograde orbits becomes smaller for deep penetrators D <C 1, because in this 
limit the binary center of mass approaches the BH in an almost radial fashion. 

If the disruption is abrupt, the ejection and capture preference could be roughly il- 
lustrated in terms of velocity (e.g Morbidelli 2006; Agnor & Hamilton 2006). The binary 
members rotate around their center of mass, such that their own motion is half of time 
with and half of time against, the motion of the center of mass r m . The net velocity of the 
members relative to the BH is accordingly increased or reduced. Since the secondary has a 
higher rotation velocity, it has more chance that the net velocity exceeds or drops below the 
escape velocity from the BH. Then, it is preferentially ejected to infinity or captured in a 
bound orbit. However, for the full discussion of the process, we also need to take into account 
the variation in the escape velocity or the variation in the potential. The displacement of 
order a in the position of each member of the binary, at a distance of about r t from the BH, 
results in a change in gravitational energy of GMa/r^ ~ ^(M/m) 1 / 3 , this is comparable 
to the variation in the kinetic energy. As we have done, it should be easier to discuss the 
overall effect in the energy domain. In our formula, the energy change fflTl) includes both 
the variation of kinetic energy and potential energy. For prograde orbits, the kinetic and 
potential terms cooperate and the net energy change is larger, the member on the "outside 
track" is expected to be ejected and its partner is captured (the "outside track" could be well 
defined especially when the orbital energy is large because the duration of the encounter is 
much shorter than the binary rotation time-scale). On the other hand, for retrograde orbits 
the variation in the gravitational energy would counteract that in the kinetic energy. 



6. Numerical Results 

In this paper we focus on results for circular coplanar binaries, though our formulae can 
be used to study the evolution of a binary with arbitrary orbital parameters. The orbit of 
a binary is assumed to be initially circular in the comoving frame of the binary center of 
mass. The center of mass of the binary is in a prograde or retrograde orbit around the BH 
(see SKR for the details of the numerical setup). 

For M/m ^> 1 the problem can be reduced to the motion of a single particle in a time- 
dependent potential ( "the restricted three-body approximation" ) described by the equations 
(JSI)- (ED and (fTUl) . The energy change, eq f[T$j) . depends only on the penetration factor D, 
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the eccentricity e, the initial binary phase 0, and the binary rotation direction. As we 
have shown, when a binary member is captured by the BH and the other is ejected, we 
can further reduce the number of the parameters by assuming e = 1 to approximate the 
additional energy. The effect of the eccentricity e 7^ 1 is taken into account through the 
non-zero initial orbital energy of the center of mass. This method (1221) will be called "the 
parabolic approximation". For a large orbital energy E (M / m) 1 ^ 3 (m/mi) (712/1712), "the 
sudden disruption approximation" (12^1) would become valid, but the three-body exchange 
reaction does not take place in this regime. In all the numerical codes the time evolution of 
objects are evaluated by using a fourth-order Runge-Kutta integration scheme. 

6.1. Energy change and probability distributions 

We have tested the restricted three-body approximation against the full three-body 
simulations of a binary evolving around a massive object (SKR). The full three-body orbit is 
accurately reproduced by the approximation equations, and the energy change, for example, 
differs at a 0.1% level when the exchange reaction happens. The comparison of the energy 
change between the full three-body and restricted three-body results is shown in figure [TJ 
Since they are in an excellent agreement, in the following discussion, we will use the restricted 
three-body approximation to test the parabolic and the sudden disruption approximation. 

In figure [U one could notice the symmetry AE(<p + ir) = —AE(<f)). The top panel shows 
that the energy change for \e — 1\/D = 0.1 is very similar to the parabolic case with the 
same D especially if we take into account the phase shift^E When the three-body exchange 
reaction happens, the value |e — 1\/D £ {mjM) x l 3 {m var jm) <C 1 should be very small, 
in such a case non-parabolic results tend toward a perfect overlapping with the parabolic 
results. In the bottom panel, the energy change in the high energy regime: (e — l)/D = 100 is 
shown. The sudden disruption approximation well reproduces the three-body results, except 
for some notable spikes. 



3 The initial distance of the binary center of mass to the BH is assumed to be r$ = 15r t for the parabolic 
calculations in the top panel of figure [TJ As long as a simulation starts at a large enough radius ro ^ r t , the 
results are largely independent of it. However, a problem arises when we compare the phase dependence. If 
we assume the same initial radius for the non-parabolic cases, it takes a slightly different time for the binary 
to reach the vicinity of the BH, the binary interacts with the BH with a slightly different binary phase. We 
have adjusted the initial radii as the periapsis passage happens at the same time (ro ~ 12.2r t for e = 0.9 and 
~ 17.4r t for e = 1.1). Since this adjustment has been done neglecting the BH tidal field, the actual phase 
at the periapsis is still different from the parabolic case. This induces the phase shifts in the figure. In the 
bottom panel of figure [U r = 15r t is assumed for the restricted three-body calculations, and in each case 
the binary phase is adjusted as the results take the maximum value at <j) = —tt/2. 
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Most binaries are disrupted at the tidal encounter with a BH. However, there are always 
finite phase regions where binaries survive (SKR). The narrow gaps in the top panel of fig[H 
correspond to such regions. Although it is not evident in the bottom panel, very narrow gaps 
exist just in the middle of the spikes. The high energy regime is usually realized with a deep 
penetrator Del. In such a case, we can ignore the self-gravity term r/r 3 in eqs ([51)-([7]), and 
free solutions are obtained. Actually, one of the free solutions: (x,y) oc (— sin/, e + cos/), 
which corresponds to the case that binary members have the same trajectory but are slightly 
separated in time, dominates around the periapsis passage (SKR). Then, the binary is once 
disrupted at the tidal encounter, but after the periapsis passage at t — 0, they come close to 
each other. If we fine-tune the initial binary phase, they form a binary again. This produces 
the narrow gaps at the spikes. If the binary phase is slightly different from the fine-tuned 
values, they almost form a binary, but they eventually break up. The additional work at 
t > due to the mutual forces between the binary members produces the notable spikes. 

The top panel of figure [2] shows the phase-averaged absolute value of the energy change 
(|AJ5|V as a function of the penetration factor D. The average is taken over the phase space 
where binaries are disrupted. We alternatively fix the orbital eccentricity or the orbital 
energy E oc (e — l)/D. For a given eccentricity, a smaller D corresponds to a larger orbital 
energy, while for a given orbital energy, it corresponds to e — > 1. For (e — 1) / D £ 1, the 
energy change remains of the order of unity as we expect from eq (123]) (see the black and 
green lines and the green crosses for the whole disruption range, and the red lines and red 
crosses for D ^ 10~ 2 ). For (e — 1)/D ^ 1, instead, (\AE\^ increases towards smaller D (the 
red lines and red crosses for D 10~ 2 ). Correspondingly, the disruption probability (bottom 
panel) becomes almost 100% for (e — 1)/D ^1. On the other hand, for the fixed energies 
(e — 1)/D = and 0.1, the disruption chance is about 80% even in the deep penetration 
limit D ^ 1 (the black and green lines). Note that in the sudden disruption approximation 
the disruption probability is 100% by definition, and we do not show it in the bottom panel. 
For our choice of (e — l)/D = 0.1, the behavior of ^A-E^ always remains very similar to the 
parabolic orbit case, even when we use the sudden disruption approximation especially for 
small D (the green crosses in the top panel). The same applies to the disruption probability 
function, but the sudden disruption approximation overestimates the disruption probability. 
As it is well known, the retrograde binaries tend to be more stable against the tidal encounter. 

For a given initial energy of a binary member / = (mi/m)E (in the following discussion 
we drop the subscript i of the initial energy for simplicity), considering the symmetry of 
AE((p), the probability P e j e (I) that the member is ejected after the tidal encounter is deter- 
mined by the fraction of the binary phase region [0, 2tt] that satisfies AE(<f)) < I, while the 
capture probability P cap {I) is determined by the fraction satisfying AE(cj)) > I. The sum 
is equal to the disruption probability: Pdi S = P e j e + P cap - Since P e j e (—I) = P cap (I), once 
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we evaluate the ejection probability P e j e (I) for negative and zero energy, the distribution 
for positive energy P eje (J) = 2P eje (0) — P e j e (—I) and the capture probability for any energy 
Pcap{I) = Peje(-I) are also obtained. The disruption probability is P^is = 2P eje (0). 

We show in figure [3] how the ejection probability of a binary member behaves for differ- 
ent penetration factors and binary orientations. Since the orbital energy and the eccentricity 
are not independent parameters when D is fixed: e — 1 oc E, we need to assume a smaller 
eccentricity, in principle, for a larger negative orbital energy I (or larger negative E) to 
estimate the energy change AE (</>), the disruption probability and the ejection probability. 
However, the energy change and the disruption probability are not so sensitive to the eccen- 
tricity, we thus evaluate them by using parabolic orbit results (the parabolic approximation) 
when the three-body exchange reaction happens. In the figure, the parabolic approximation 
calculations (the lines) are in a good agreement with the restricted three-body calculations 
(the circles). For prograde orbits, the energy change AE(<f>) is more sensitive to the phase 
4> around zero points, only very narrow phase regions satisfy \AE\ 1 (see fig [1] top panel), 
then the ejection probabilities have a plateau around 1 = 0. The long low energy tail of the 
prograde D = 0.1 case (the green solid line) reflects the fact that the energy gain would be 
quite large for a narrow binary phase region for this case. The ejection (capture) probability 
for — oo < I < oo is a monotonically increasing (decreasing) function of the energy I. The 
probabilities rapidly change around |/| ~ 1 and their function form depends mainly on D 
and the binary rotation direction. 



6.2. Probability distributions for various semi-major axes 

When we discuss an actual astrophysical system, it is more physically intuitive to use 
the semi-major axis rather than the energy I or E. By specifying the semi-major axis and 
the mass ratios, the ejection (or capture) probability and the energy after the disruption can 
be evaluated. The orbital energy E is given by 

p _ o i ( r a/ r t \ ~* ( M/m \ 1/3 / m 1 /m \ ~ X ( m 2 /m \ -1 

E -~ 2 - 7 \im) \-w~) \~3jr) ■ [ } 

M/m = 10 6 and mxjm^ = 3 will be assumed in this section. 

Figure H] shows the characteristics of the ejection process when a binary approaches 
the BH in elliptic orbits. To evaluate the ejection probability (the top panel), we have not 
distinguished which member is ejected. Since both members are never ejected together when 
E < 0, it is just the sum of the ejection probabilities of the two members. For parabolic 
orbits, one of the members is always ejected if the binary is disrupted, then the black solid 
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line indicates the disruption probability. As r a /r t becomes smaller, the deviation of the 
orbit from the parabolic one becomes larger. The effect of the non-zero orbital energy is 
expected to become significant around r a /r t ~ ( y M/m) 1 ^ 3 (m/m 2 ) ~ 400, and the ejection 
probability rapidly decreases. The middle panel indicates which member in a binary is ejected 
more frequently. There is no preference in the parabolic case. However, the secondary is 
preferentially ejected for small r a /r t . The phase-averaged ejection energy is shown in the 
bottom panel. For parabolic orbits this is equivalent to the phase-averaged energy change 

(my 

If we consider P cap (I) = P e j e (—I), it is possible to interpret figure 0] as the capture 
probability and capture preference for the hyperbolic orbit cases with the same semi-major 
axises in absolute value. The phase-averaged energy (the bottom panel) gives the absolute 
values of the averaged capture energies. 

7. Irregular satellites around giant planets 

Over 150 satellites are orbiting the giant planets in the Solar system. About one-third 
of these are classified as regular, with nearly circular and planar orbits. The majority of 
the satellites, however, are irregular ones which are more distant from their planet and 
typically have larger eccentricities and/or inclination. Interestingly, a large fraction of the 
irregular satellites orbit their planet in the retrograde direction (Jewitt & Haghighipour 
2007). Triton, Neptune's largest moon, is among them. Because of the retrograde orbit and 
composition similar to Pluto's, Triton is thought to have been captured from the Kuiper 
belt. Recently Agnor and Hamilton (2006) demonstrated that the gravitational encounter 
between Neptune and a binary system, which had included Triton as a member, is an effective 
mechanism to capture Triton. Since the mass ratio between Neptune tjin ~ 17.2M e and 
Triton uit ~ 3.58 x 10~ 3 M e is reasonably large m^lvnT ~ 4800, we here revisit the capture 
process from the point of view of the restricted three-body problem. 

We consider an encounter between Neptune and a binary system with m\ = m-i = 
O.lrriT and a = 20Rt- The binary is assumed to be in a prograde hyperbolic orbit with the 
periapsis r p = 8Rjy where Rt ~ 1.35 x 10 3 km and Rn ~ 2.46 x 10 4 km are the radii of 
Triton and Neptune, respectively. This set of the parameters are identical to what Agnor 
and Hamilton (2006) have assumed to obtain their figure 2. Since the tidal radius is about 
r t ~ 18Rn, the penetration factor for the orbit is D ~ 0.45. For a given encounter velocity 
at infinity Voo, we evaluate the binary disruption chance and the capture probability of the 
binary members. The velocity at infinity Vqo — 1 km/s is related to the semi-major axis of 
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the hyperbolic orbit as 

K/rtl = 16 (iTrnTi) G^fr) (li) (i^w) ' (26) 

which is about the critical value \r a /r t \ ^ (M/m) 1 / 3 ~ 16 for the three-body exchange re- 
action. The numerical results are shown in the figure \5[ and the relevant case (the solid 
lines in the top panel) is in good agreement with the figure 2 in Agnor and Hamilton (2006). 
Since the disruption chance depends mainly on D, for the fixed D, it is constant ~ 97%. 
For high velocities, the probability for the primary drops sharply, while it is still ~ 50% for 
the secondary because of its smaller initial energy /. If we assume a larger secondary mass 
m 2 = rriT/2, the energy change AE at the encounter would be larger. The 50% capture 
probability of the primary now extends to higher velocities, while the capture rate of the 
secondary is similar because the primary mass is fixed (the dashed lines) . The black dashed- 
dotted lines indicate the capture probabilities obtained by using the parabolic approximation 
for m 2 = O-lrriT- Even at the high velocity Voo ~ 1.6 km/s or equivalently high eccentricity 
e — 1 = D(m/M) 2 ^ 3 (avl /Gm) ~ 7 x 10~ 2 , it reasonably agrees with the restricted three- 
body results. In the bottom panel, the capture probability is shown as a function of a scaled 
velocity (m/mT) 1 ^ 6 (m par /mT)~ 1 ^ 2 v 00 = (m^ / 'm^) 1 ^ \2IGrriT / 'a) 1 / 2 . Although the difference 
among them is small, if the parabolic approximation is used for all the cases, the solid and 
dashed lines should perfectly overlap with the black dashed-dotted line. 

At the disruption, each member is captured with 50% probability if the initial energy 
I is close to zero. Using the relation P cap {I) = P e je{—I), figure [3] also shows that at high 
energies the capture is very rare. A transition occurs at an intermediate energy as shown in 
figure El For D = 0.45 prograde orbits, the capture probability sharply drops from ~ 50% 
around I ~ 1.4. The critical value of at which the capture probability of a binary member 
mi drops from ~ 50% is given by 

l2IGm 2 [M\ 1/& t fl\ 1/2 f m 2 \ 1/2 / a \" V V M/m . 

(27) 

This estimate well explains the critical velocities in figure |5j If Triton has a heavier compan- 
ion and/or a member of a harder binary, the critical velocity could be higher, provided the 
binary is disrupted by Neptune: D ^ 1. The semi- major axis of the captured member mi is 

2\E 1 \\m 2 ) \MJ Nl 11 \20R T J \4.80 x 10 4 J \4.37 x 10 3 y ; 

where \E\\ is about \AE\ ~ 1. After capture, the orbit of Triton needs to shrink to the 
current observed r a ~ 14R^ and e ~ 10 -5 either through tides or other means (Correia 
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2009; Nogueira et al. 2011). The critical velocity and the semi-major axis for the secondary 
m 2 are obtained by exchanging the subscript 1 and 2 in eqs (1271) and (128]) . 

Although a detailed modeling of the Solar system is beyond the scope of this paper, 
the Sun can not be generally ignored in the capture process (e.g. Philpott et al. 2010 
Gaspar et al. 2011). Since the tidal radius r t ~ 18i?Ar(m/m;r) _1 (a/20-RT) is much smaller 
than the radius of Neptune's Hill sphere r# ~ 4700-Rat, Neptune completely dominates the 
attraction of a binary during the tidal encounter. Tidal effects from the Sun are negligible for 
the disruption process itself, which is the main focus of the paper. However, the permanent 
capture of Triton requires an additional condition that after the tidal breakup a capture orbit 
should be well inside the Hill sphere. Detailed numerical simulations show that retrograde 
satellites (those orbiting in the opposite sense as Neptune orbits the Sun 0) of Neptune are 
more stable than prograde ones, and they are stable to distance of r stab i e ~ 0Ar H (Nesvorny 
et al. 2003; Holman et al. 2004). Capture orbits are more eccentric than assumed in the 
numerical simulations, the stability region might be slightly smaller. Considering that the 
apocenter distance is about twice the semi-major axis ( )28|) for highly eccentric orbits, the 
binary which had included Triton as a member should satisfy a relation, 

/ a \ / m 2 \ 1 / m \ ^ 3 < / r sta bie \ 
\20^) V0.1m T y \l.lm T ) ~ ' \0.3r H ) ' 

where we have used \AE\ ~ 1. A smaller value of the ratio a/rri2 gives a smaller apocenter 
distance and a higher critical velocity (|27l) . As shown in the top panel of figure [2] (the black 
solid and dashed lines), the phase-averaged value (|Ai?|) is slightly larger than unity for 
prograde binaries, and slightly smaller for retrograde binaries when the disruption probability 
is high. For a given encounter velocity, a retrograde orbit of the binary center of mass around 
Neptune with a prograde binary rotation is optimal to produce a stable capture orbit. 

The ellipticity ~ 1.7 x 10~ 2 of Neptune induces a small deviation in its own grav- 
itational potential from the point-mass estimate especially in a non-planar configuration. 
However, the deviation in the attraction force AF/F and the tidal acceleration Aa t /at is of 
order of (Rn /rt) 2 £N ~ 5 x 10 -5 at the tidal radius. The effect is clearly negligible. 



In the three-body encounter discussion, prograde motion means that the binary center of mass is orbiting 
around a massive object in the same sense as the binary members rotate around their center of mass. Then, 
a moon can be captured in a retrograde orbit around Neptune (the orbit is in the direction opposite to the 
rotation of Neptune) after the tidal breakup of a prograde binary (the angular momentum of the binary 
around Neptune and of a binary member around the binary center of mass are aligned). 
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8. Conclusions 

We have discussed how the members of a binary are ejected or captured after a tidal 
encounter with a massive object. We have shown that the ejection and capture dynamics can 
be well understood in the framework of the restricted three-body approximation. When the 
three-body exchange reaction happens, the orbit of the center of mass of the binary should 
be almost parabolic 1 1 — ej ^ D(m/M) 1 ^ 3 (m par /m) <C 1 or equivalently the semi-major axis 
is large \r a /r t \ (M /m) l l 3 {m/m par ). The essential quantity to characterize the disruption 
process is the energy change AE at the encounter, which practically depends only on three 
parameters: the penetration factor D = r p /r t , the binary rotation direction, and the binary 
phase 0. Except for the phase shift of n, the energy change is exactly the same for the two 
members with arbitrary mass ratio. 

In principle, for any positive (negative) orbital energy E of the center of mass of a binary, 
the heavier member has more chance to be ejected (captured), because it carries a larger 
fraction of the orbital energy. However, if the orbital energy is close to zero, the difference 
between their ejection (capture) probabilities becomes small, and there is practically no 
ejection and capture preference. For a parabolic orbit, each member is ejected in exactly 
50% of the cases. The preference becomes significant when the absolute value of the energy 
\E\ is comparable to the typical energy change (Gm 1 m 2 / 'a)(M / 'm) 1 / 3 . On the other hand, 
if \E\ is much larger than the typical energy change, they are both ejected for E > or 
captured for E < 0, and there is no ejection or capture preference. 

Corresponding to the typical energy change, we can define a critical encounter velocity 

^oo, crit = 

(2IGm par / 'a) l / 2 {M /m) 1 / 6 for the capture process where the critical initial energy 
/ ~ 0.5 — 1.5 mainly depends on the penetration factor and the orientation of the binary. 
Since the distribution of the energy change is the same for the two members, the secondary 
star has a higher critical velocity. The tidal encounter might be responsible for the capture 
of not only Triton but also other irregular satellites in the Solar system (Morbidelli 2006, 
but see also Vokrouhlicky et al. 2008; Philpott et al. 2010. The stability of capture orbits 
for the solar perturbations also should be taken into account). However, the other irregular 

1/2 

satellites are much lighter than Triton. The capture mechanism w 00j crit oc m par requires that 
these irregular satellites have been in a binary with a very massive partner and that the less 
massive member has been predominantly captured. 

If a binary has a large orbital energy E 3> (M /m) 2 ^ 3 (Gm 2 / a) , the disruption would 
happen slightly inside the usual tidal radius at which the BH tidal forces become comparable 
to the binary mutual gravity forces. The abrupt disruption approximation provides a good 
estimate of the energy change, and the typical energy change in the high energy regime is 
much larger than in the low energy regime. However, the change is not significant compared 
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to the initial energy E, the tidal encounter is not an efficient acceleration process anymore. 
In the high energy regime, both members are ejected after the BH encounter. 

In recent years, observations have identified a remarkable number of hypervelocity stars 
(Brown et al. 2009; Tillich et al. 2011). The distribution of line-of-sight velocities of 
hypervelocity star candidates shows a long tail in the high velocity region (v > 275 km/s) 
which includes comparable numbers of unbound and bound stars assuming that the escape 
velocity of the Galaxy at 50kpc is ~ 360 km/s (Brown et al. 2009; Kenyon et al. 2008). 
Since the initial orbital energy is negligible when the three-body exchange reaction happens, 
the velocity distribution should reflect the distribution of the energy change AE at the tidal 
encounter (and the Galactic potential) with its dependence on the periapsis distance, the 
binary orientation and phase. 

Binaries are supplied to the BH at the Galactic center predominantly from its radius 
of influence ~ a few pc or even beyond it (e.g. Perets et al. 2007; another interesting 
possibility is that they might arise from the nuclear stellar disk. e.g. Madigan et al. 2011). 
Those that come from about the radius of influence are on elliptical orbits. Since the radius 
of influence is much larger than the tidal radius, the preferential ejection for high and low 
mass members (m x > m 2 ) is therefore irrelevant for stellar binaries (mi/m 2 ~ 10) if the 
binary semi-major axis is smaller than ~1AU. However, stars with planets imxjm^ 10 3 ) 
of semi-major axis of 1AU will predominantly eject the planets. 

We thank the referee, Hagai Perets, for his constructive comments, and Mostafa Ahmadi, 
Phil Armitage, Witold Maciejewski, and Ruth Murray Clay for useful discussion. This 
research was supported by STFC, ERC, and IRG grants, and Packard, Guggenheim and 
Radcliffe fellowships. 
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(J)/tc 

Fig. 1.— Energy change as a function of <fi: top panel (£) = 1): the restricted three- 
body approximation for e = 0.9 (green solid line), 1 (red solid line), and 1.1 (black solid 
line). The full three-body calculations (circles). Bottom panel (D = 10~ 3 and e = 1.1): the 
restricted three-body approximation (black solid line), the full three-body calculations (black 
circles) and the sudden disruption approximation (red solid line). In the full three-body 
calculations, M/m = 10 6 and mi/mj = 3 are assumed and the energy change is evaluated as 
AE = (m\/m)Ei — ^/m)^. Prograde orbits are assumed for all the calculations. Energy 
is in units of {Gmxm^j a)(M /m) 1 ^ 3 . 
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Fig. 2. — Phased- averaged energy change (top panel) and disruption probability (bottom 
panel) as a function of D. Results for the restricted three-body approximation are shown 
for e = 1 (black lines), e = 1.01 (red lines) and (e — 1)/D = 0.1 (green lines). Prograde 
(solid lines) and retrograde (dashed lines). The sudden disruption approximation is shown 
for prograde orbits with e = 1.01 (red crosses) and (e — 1)/D = 0.1 (green crosses). Energy 
is in units of (Gmim2/a)(M/m) 1 ' 3 . 
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Fig. 3. — Ejection probability as a function of / = (mi/m)E. The parabolic approximation 
is shown for D = 1 (black line), 10 -1 (green lines) and 10 -2 (red lines). The solid and dashed 
lines indicate prograde orbit and retrograde orbit results, respectively. The circles show the 
restricted three body results for the corresponding cases with M/m = 10 6 and rrti/m = 1/4. 
Energy is in units of (Gmiiii2/ a)(M /m) 1 / 3 . 
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Fig. 4. — Ejection probability, preference and energy. Parabolic orbit (black solid) , \r a /r t \ = 
500 (red solid), 200 (green solid), 130 (red dashed) and 100 (green dashed). The parabolic 
approximation is used. The ejection preference is the ratio of the ejection probabilities (the 
primary star/the secondary star). The ejection energy is evaluated by taking the phase 
average of all the ejected cases. Energy is in units of (Gmim 2 /a)(M/m) 1 / 3 . 
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Fig. 5. — Top panel: The restricted three-body approximation: binary disruption chance 
(black solid), capture chance for m 2 = O.lm^: the primary (red solid) and the secondary 
(green solid), capture chance for m 2 = m^/2: the primary (red dashed) and the secondary 
(green dashed). The parabolic approximation for m 2 = O.lmr (black dashed-dotted). Bot- 
tom panel: Capture chance as a function of the scaled velocity. D ~ 0.45 is assumed. 



